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Abstract 
Since  Pakistan  has  varying  climates  and  terrains,  the  Pakistan  Army  ro- 
tates its  units  between  peacetime  locations  so  that  no  unit  endures  inequitable 
hardship  or  enjoys  unfair  advantage.  Arm]/  policy  spec/In  s  strict  constraints 
on  mnt  rotations,  such  as  tin  lenqlh  of  a  unit's  stay  in  any  location,  tin 
number  of  units  moving  at  any  tune,  and  the  allowable  replace  menls  for  any 
moving  unit.  Scheduling  rotations  manually  in  accordance  with  these  rules, 
as  is  currently  practiced,  is  extremely  difficult  and  lime-consuming.  This  pa- 
per presents  an  integer  programming  model  thai  finds  feasible,  minimum-cost 
schedules  for  the  Pakistan  Army's  desired  planning  horizons.  Ihe  model  also 
ensures  thai  the  units  are  positioned  at  Ihe  end  of  the  planning  horizon  so 
that  feasible  schedules  exist  for  future  planners.  The  model  is  implemented 
with  commercially  available  optimization  software.  Schedules  an  obtained 
feir  realistic  lest  problems  in  less  than  an  hour  on  a  personal  computer. 


The  Pakistan  Army  peacetime  rotation  problem,  as  described  in  the  ab- 
stract, was  chosen  by  Colonel  Baig  as  a  research  topic  for  a  masters  thesis. 
I3aig  [1992],  directed  by  his  co-authors.  lie  received  guidance  on  problem  def- 
inition from  Brigadier  General  Maroon  Basin r  Khan,  who  was  then  serving  as 
the  Army  Attache  at  the  Pakistan  Kmbassy  in  Washington.  The  model  pie 
sented  in  this  paper  captures  all  aspects  of  the  peacetime  rotation  problem 
as  described  by  the  general. 

The  peacetime  rotation  problem  is  a  unique  scheduling/timetabling 
problem.    It.  shares  the  diversity  and  large  number  of  constraints  typically 
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found  in  the  construction  of  sport  schedules.  (E.g.,  sec  Andruc  and  Coromi- 
nas  [1989]  for  the  1992  Summer  Olympic  Games,  Shell  [1985]  for  the  Na- 
tional Basketball  Association,  Cain  [1977]  and  Shell  [19X5]  for  rVlajor  League 
Baseball,  and  Ferland  and  Pleurent  [1991]  for  the  National  Hockey  League). 
However,  the  peacetime  rotation  problem  has  relatively  few  events  compared 
to  the  hundreds  or  thousands  of  games  that  make  up  a  typical  spoils  sched- 
ule. An  exact  procedure  is  therefore  undertaken  in  this  paper  instead  of 
the  heuristic  procedures  used  successfully  for  large  scale  timetabling  prob- 
lems Carter  [1986].  The  following  sections  present:  I)  policies  defining  the 
problem,  2)  the  integer  programming  formulation,  3)  model  refinements  for 
increased  tract  ability,  1)  computational  experience,  and  5)  conclusions. 

1      Pakistan  Army  Peacetime  Rotation  Policies 

The  Pakistan  Army  classifies  inilil  ary  local  ions  into  Peace  Areas  (  PAs).  Semi- 
Hard  Areas  (SIIAs)  and  Hard  Areas  (HAs).  This  classification  accounts 
for  Pakistan's  diverse  terrain,  ranging  from  desert  to  lofty  mountains,  and 
temperatures,  ranging  from  below  freezing  to  above  10  degrees  Celsius.  The 
geographic  classifications  also  account  for  available  facilities  al  the  locations 
and  the  proximity  to  major  metropolitan  areas.  To  ensure  thai  personnel 
serve  equally  in  all  three  areas,  unit  personnel  are  rotated  between  locations 
during  peacetime  on  a  regular  basis. 

Military  units  in  the  Pakistan  Army  are  classified  in  three  ways:  aceoi  ding 
to  their  operational  role  (strike  oi  defensive),  according  to  their  functional 
role  (Armor,  Artillery,  Engineers,  Infantry,  Signals,  Supply,  etc),  and  ac- 
cording to  whether  their  equipment  is  supplied  by  Eastern  or  Western  bloc 
nations.  (See  Cohen  [19X1]  for  the  roles  of  the  East  and  West  in  equipping  the 
Pakistan  Arm}-.)  Separating  the  units  into  distinct  categories  based  on  these 
three  attributes  simplifies  subsequent  mathematical  analysis  and  is  useful  for 
elucidating  rotation  policy. 


The  General  Headquarters  oi  (lie  Pakistan  Army  lias  established  rotation 
policies,  whoso  salient  features  are  as  follows. 

1.  A  unit  ran  movo  only  il  replaced  by  a  unit  ol  the  same  category,  i.e., 
having  the  same  operational  and  functional  roles,  and  operating  the 
same  type  of  equipment.  This  results  in  a  natural  division  of  the  overall 
problem  into  separate  rotation  problems  for  each  category.  The  largest 
category,  currently,  has  <ST  units. 

2.  If  a  unit  at  location  A  moves  to  location  13,  then  a  unit  of  the  same 
category  at.  location  H  must  move  to  location  A  in  the  same  year. 
This  policy,  referred  to  as  mutual  replacement,  simplifies  transfer  ol 
operational  and  administrative  responsibilities  at  both  locations. 

3.  Equipment,  does  not  move,  only  personnel  do. 

4.  No  more  than  one  unit  can  move  from  the  same  brigade  in  the  same 
year.  A  brigade'  is  composed  of  thiee  units  ol  the  same  category. 

5.  Some  units  (such  as  Engineers,  Signals,  Reconnaissance  and  Support 
Battalions)  do  not  belong  to  a  brigade.  They  fall  under  direct  control 
of  a  Division.  'I  heir  can  be  no  more  than  one  of  these  units  moving 
from  the  same  location  in  the  same  year. 

6.  Each  unit's  tenure  requirement  varies  by  location  classification  as  fol- 
lows: 

•  Peace  Area  (PA)  5-7  years, 

•  Semi- Hard  Area  (SUA)  2-4  years, 

•  II ard  Area  ( 1 1 A )  1  - 3  years . 

7.  An  individual  unit  must  rotate  according  to  the  cycle  ol  locations: 
PA  -»  SUA  -»  PA  -»  II A  ->  PA.  as  shown  in  Figure  I. 


Figure  1:  Each  unit  must  rotate  through  Peace  Areas.  Semi-Hard  Areas  and 
Hard  Areas  in  the  indicated  order.  Alternate  Peace  Areas  in  the  cycle  must 
be  different.    The  length  of  stag  in  each  location  mast  Jail  within  a  prescribed 

range. 


S.  There  is  no  restriction  placed  on  which  SHA  or  HA  locations  a  unit 
visits  on  its  cycle.  However,  a  unit  must,  not  return  to  its  previous  PA 
location. 

Currently,  peacetime  rotation  schedules  are  developed  manually  by  plan- 
ners at  the  General  Headquarters  of  the  Pakistan  Army.  A  live-year  schedule 
is  developed  on  a  yearly  basis.  Due  to  the  large  problem  size  and  complex; 
policy  structure,  the  current,  method  suffers  from  the  following  drawbacks: 

•  The  units  may  be  positioned  at  the  end  of  the  planning  horizon  in  a  way 
that  precludes  future  schedules  from  satisfying  the  policy  constraints. 

•  The  schedule  requires  hundreds  of  man-hours  to  develop. 

•  It  is  difficult  to  evaluate  proposed  policy  changes. 

•  The  schedule  may  not  be  developed  impartially. 


•  The  schedule  may  incur  excess  cos!  by  transferring  units  more  oil  en  or 
over  greater  distances  than  necessary. 

These  limitations  motivated  thr  development  of  an  integer  programming 
model  to  assist  with  rotation  scheduling. 

2      Peacetime  Rotation  Model 

Our  integer  programming  model  for  the  Pakistan  Army's  peacetime  rotation 
problem  minimizes  the  total  cost  associated  with  all  scheduled  moves,  while 
ensuring  that  all  policy  constraints  are  satisfied.  The  model  is  valid  for 
time  horizons  up  to  15  years,  which  satisfies  the  Pakistan  Army's  planning 
requirement  for  five-year  rotation  schedules. 

•  Indices: 

i,i'  units, 

/,  /'  locations, 

/,/'  years. 

•  Given  Data: 

null,  the  initial  location  for  unit  i. 

mini  minimum  stay  allowed  at  location  /. 

max i  maximum  stay  allowed  at  location  /. 

mcostmi  movement  cost   for  unit  i  from  /  to  /', 

stayn  number  of  years  unit  i  lias  been  at   location  /  at   the  start  ol  the 

first  time  period. 

V         set  of  peace  area  locations. 

H  set  of  hard  area  locations, 


61         set.  of  semi-hard  area  locations. 

•  Derived  Sols:  The  following  sets  arc  used  to  enforce  rotation  policies. 
Their  derivation  is  discussed  in  section  3. 

Uii't  set  of  all  units  eligible  to  move  from  location  /  to  /'  in  year  t, 

Fm         set  of  all  possible  locations  from  which  unit  i  could  have  moved  il 
it  arrives  at  local  ion  /  in  year  /, 

I'm  set-  of  all  possible  local  ions  to  which  unit  ?  can  move  il  sit  uai.ed  at 

location  /  in  year  /. 

•  Decision  Variables: 

T,ii't  ]    if  ti ii il    ?  moves  from  location   /  to  location  /'  in  year  /,  and  0 

otherwise. 

•  Formulation: 


minimize  E  E  E  E  nicest 


i',ii'riii'i 


r         l        l'        1 


mn.ri  —uta ;/,; 

E      E  •'■'"''  =  '  V//  ■"■'■  '  =  ,mt,< 

tz=mini—.ilny,i  /'£  I ",,, 

E  ■'»'»  =  E  mi'*  v//,/ 

iet',,|,  feu,,,, 

E  E  *»i't  <  1  v// 

t  H  vuiTi 

E  E  '.■/'//'<    E     E  *.■//'«•  V/// 

t  —  mitif 

E  •r'»''  ^    E     E  ■r'''"'  V/// 


(I) 


(3) 


(5) 
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££I>.-'/'f  <  "  v;  (r>) 

££5>f7/i,  <  I     V;  (7) 


L  E  ■'•-"''  <  i  v/ 


z. 


EE^"  +  EE^™<  j  v?./g  r  (in) 

few    '  fes    ' 

EEEE''>'  =  EEEE^'  (|]) 

( 'onstra.int  ( 1 )  ensures  carl i  unit  s  first  move  is  completed  within  (lie  min- 
imum and  maximum  tenure  requirements  at  its  initial  location.  Constraint 
(2)  enforces  the  mutual  replacement  policy.  Constraint  (•'?)  allows  no  more 
than  one  unit  to  move  from  (lie  same  location  in  the  same  year.  II  two 
or  more  brigades  have  the  same  geographical  location,  this  location  is  bro- 
ken into  separate  brigade  locations.  Thus,  constraint  (3)  also  enforces  the 
restriction  that  no  mote  than  one  unit  can  move  from  the  same  brigade. 

Const  taints  (  1 )  and  (.r))  control  the  tenure  requirement  s  for  all  moves  fak- 
ing place  after  the  first  move.  Constraint  (I)  restricts  the  number  ol  times 
a  unit  enters  location  /  up  to  year  /  to  be  less  than  or  equal  to  the  number 
of  times  it  leaves  th*1  location  up  to  year  /  +  ina.r\.  These  constraints  ensure 
that  a  unit  leaves  a  location  in  year  /  only  il  it  arrived  at  that  location  in  an 
appropriate  earlier  year.  (They  can  alternatively  be  formulated  with  noncu- 
mulat  ivo  inequalities,  but  the  cumulat  ive  form  provides  better  computational 
performance.)  Constraint  (-r>)  requires  a  unit  leaving  location  /  in  yeai  /  to 
have  arrived  there  /  —  maxi  to  /  —  m/vi  years  earlier.    Constraints  (  I)  and 


(!">)  cannot,  prevent  a  unit  from  moving  away  from  I, lie  same  location  more 
Mian  once  in  I. lie  last  few  years.  Constraints  (0)  -  (0)  eliminate  this  problem 
by  allowing  each  unit  at  most  one  move  between  difTcrcnl  areas  (hiring  the 
planning  horizon  (a  valid  restriction  for  planning  horizons  of  15  years  or  less). 

Constraint  (10)  prevents  a  unit  from  moving  to  the  same  PA  more  than 
once.  If  the  planning  horizon  is  less  than  seven  years,  this  constrain!  is 
unnecessary. 

Constraint  (11)  states  the  total  number  of  moves  from  II  As  to  PAs  should 
be  equal  to  the  total  number  of  moves  from  SIIAs  to  PAs.  This  constraint 
helps  position  units  appropriately  at  the  end  of  the  planning  horizon,  as 
discussed  further  in  section  '.]. 

3      Model  Refinement 

Even  though  we  can  treat  the  Pakistan  Army's  rotation  scheduling  problem 
with  separate  models  for  each  operational/functional  area,  the  independent 
models  can  still  be  quite  large.  The  largest  operational/functional  area  has 
87  units  spread  over  30  locations.  Straightforward  application  of  the  pre- 
ceding model  for  these  units  over  an  eight-year  horizon  would  require  over 
000,000  binary  variables.  Fortunately,  characteristics  ol  the  policies  can  be 
exploited  to  identify  many  impossible  unit  movements  and  eliminate  the  cor- 
responding variables.  The  transition  eligibility  parameter  OKiWi  *s  defined 
for  this  purpose.  It  has  value  1  if  and  only  if  unit  i  is  eligible  to  move  from 
location  /  to  /'  in  year  /.  The  idea  is  to  make  this  parameter  zero  as  often  as 
possible  without  sacrificing  model  fidelity  or  optimality.  The  sets  Uwi,  J',n. 
Fni  defined  in  section  2  are  readily  generated  from  this  parameter.  The  logic 
for  deriving  OK,wt  is  as  follows. 

1.  Rotation  policies  and  tenure  limits  severely  restrict  the  set  of  allowable 
moves.  For  example,  a  unit  that  has  moved  from  a  HA  and  has  been 
at  a  PA  for  three  years  can  not  possibly  move  for  two  more  years  and 
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then  only  to  a  SUA.  Subsequent  moves  are  similarly  restricted.  I  hose 
situations  are  detailed  below  in  four  exhaustive  cases.  Let  OP  A,  denote 
the  current-  or  most  recent  FA  location  for  unit  /  in  the  initial  year, 
Stay  —  m\n{$tayn,mirii}  for  the  unit's  initial  location  /,  and  let  "\" 
denote  set,  exclusion. 

a )  If  unit  i  is  current  I  v  located  at  a  PA  and  previously  served  in  a  SI  1  A, 


OK, 


ii't 


I    if  I  G  P. 

i  if  i  g  n, 

1    if  I  e  {P\OPA, 

0    ill  Ik  i  in  si  . 


I'  en, 

I'  g  {p\opa, 

V  g  5, 


6  <  /  +  Slay  <  8. 

7  <  /  +  Stay  <  11, 
12  <  /   I-  Slay  <  IS. 


b)  If  unit  i  is  currently  located  fit.  a  PA  and  previously  served  in  fi  II  A. 


OK  ,ii' i  =  I 


I    /./'    /  G  P. 

I    if   I  G  S, 

I    if   if  {P\()PA, 

0   otherwise . 


/'  G  5, 

/'  G  \P\OPA, 

r  g  k 


(i  <  /  +  Slay  <  8. 
8  <  /  +  .S'/r/iy  <  12. 
13  <  /  +  Stay  <  19, 


c)  If  unit  ?  is  currently  located  at  a  SUA. 


OKuvt  =  { 


I    if  /G  5,  /'G  {P\0/M, 

1    ?/  /G  {P\OPA,}.    V  e  V. 
\   if  leH,  /'G  P. 


I)     r>// 


lerinsi  . 


:{  <  /  +  5/fl.y  <  5, 
8  <  /-}  57-7?/  <  12, 
<)  <  /  -|-  Slay  <  15, 


(I)  II  unit  7  is  currently  located  ri  1  a  IIA. 


OKuvi  =  < 


I    if  I  £H,  V  £{V\OPAi},  2<t  +  Slay< 

I    if  I  6  {V  \  OP  A,},  V  G  5,  7  <  H  tto?/  < 

I    ?./'   /  e  5,  /'  £  'P,  9  <  /  t  Stay  < 
0   otherwise. 


1. 
I  I. 
15. 


For  a  detailed  example,  suppose  unil  i  is  currently  located  at  a  PA  and 
its  old  location  is  a  IIA  (case  b).  This  unil  can  move  to  a  SUA  after 
completing  5  to  7  years  of  stay  at  the  PA,  so  it  is  eligible  to  move  only 
when  /  -f  Stay  is  6,  7  or  8  years.  The  unit's  next  move,  to  any  PA 
except  OP  A,  must  lake  place  2  to  I  years  later,  i.e.,  in  one  of  the  years 
8  through  12. 

2.  The  mutual  rotation  policy  also  helps  eliminate  many  variables.  Con- 
sider a  unit  at  location  /  that  is  eligible  to  move  to  I'  in  year  /  after 
executing  step  1 .  That  move  ran  be  scheduled  only  il  another  unit  is  el- 
igible to  move  from  local  ion  /'  to  local. ion  /  in  the  same  year.  Therefore, 
any  OI\ju>,  thai  was  1  after  step  1  is  changed  to  0  unless: 

£OAVn,  >  I- 

3.  We  can  extend  the  idea  of  step  2  to  subsequent  moves,  for  unit  >  to 
be  eligible  lor  a  move  to  location  /'  in  year  /,  there  must  be  another 
unit  ?'  eligible  to  replace  unit  ?  at  location  /'  between  /  \  minii  and 
/  -f  max\>.  Therefore,  any  OKmn  thai  remains  I  after  the  first  two 
steps  is  changed  to  0  unless: 

EI    E    o/w>i. 
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YearO 


Year  3 


Year  6 
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Figure  2:  A  feasible  six-unit  rotation  pattern  that  can  be  extended  indefinitely. 
Nodes  in  the  figure  represent,  units  undergoing  moves  in  the  indicated  year. 
Edges  represent,  mutual  replacement. 

3.1      Conditions  for  Future  Feasibility 

The  Pakistan  Army  needs  to  ensure  that  units  are  positioned  at  the  end 
of  the  planning  horizon  so  that  feasible  schedules  exist,  for  future  planners. 
Sufficient  conditions  are  developed  for  this  purpose.  These  conditions  are 
explained  with  the  help  of  Figure  2. 

Figure  2  shows  an  indefinitely  repeatable  rotation  schedule  that,  satisfies 
all  restrictions  of  the  rotation  policy  outlined  in  section  L.  This  figure  con- 
tains 4  PAs,  1  HA,  and  1  SHA  locations  where  each  location  has  one  unit. 
The  circles  contain  unit  identifiers  and  the  arcs  between  the  circles  indicate 
units  exchanging  locations.  For  the  units  initially  located  al  PAs  (1,  2.  3, 
and  6),  their  last  area  is  shown  with  H  (for  HA)  or  S  (for  SHA). 


II 


It  is  possible  to  add  units  in  multiples  of  6  (4  PA  units,  I  IIA  unit,  I  SUA 
unit)  up  to  the  maximum  of  three  at  each  location,  Kach  (i  additional  units 
(1  at  each  existing  location)  could  be  feasibly  added  to  the  above  schedule 
by  allowing  moves  in  similar  3  year  increments  starting  at  year  I  for  the  first 
G  units  and  at  year  2  for  the  second  (i.  New  locations  with  units  following 
the  same  pattern  could  also  be  added. 

The  positioning  of  units  in  Figure  2  obeys  the  following  onditions  at  all 
times. 

Condition  1  The  total  number  of  units  in  PAs  is  twice  the  number  of  units 
in  II As  and  SIIAs. 

Condition  2  The  number  of  units  in  IIAs  and  SIIAs  is  equal. 

Condition  3  Half  the  units  in  PAs  previously  served  in  SIIAs  and  half  pre- 
viously served  in  11  As. 

The  three  conditions  above  are  not  necessary  to  guarantee  the  existence 
ol  a  feasible  solution  in  the  future.  However,  as  proven  in  the  appendix,  if 
the  conditions  are  satisfied  for  the  six  years  preceding  the  current  horizon 
and  are  enforced  throughout  (he  horizon,  then  future  feasibility  is  guaranteed. 
Unfortunately,  some  operational/functional  areas  have  not  always  conformed 
to  the  conditions.  Therefore,  the  test,  problems  of  section  I  are  solved  with 
constraint  (11)  instead  of  the  more  restrictive  form  of  the  constraint 

£E£*.«"  =  X;EX>.'»"  v'  (J 2) 

which  would  guarantee  feasible  rotations  in  the  future  wit  h  appropriate  start  - 
ing  conditions.  Though  not  guaranteed,  constraint  (II)  empirically  yielded 
rotation  schedules  with  ending  conditions  that  allowed  future  feasible  rota- 
tions. 
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4      Computational  Experience 

The  Pakistan  Army  peacetime  rotat ion  model  is  implemented  in  the  (leneral 
Algebraic  Modeling  System,  (JAMS,  [1992]  and  solved  using  X  A  [1987]  and 
OSL  [1991].  All  computational  results  are  obtained  using  a  480/33  personal 
computer  with  16  megabytes  of  HAM.  The  goal  of  our  work  was  to  develop 
an  implenienta.ble  scheduling  model.  We  believe  our  choice  of  commercially 
available  software  and  a  personal  computer  represents  the  best  chance  of 
implementation  for  the  following  reasons: 

1.  The  Pakistan  Army  can  implement  the  model  for  a  reasonable  cost. 

2.  The  software  is  stable,  well  documented,  and  the  user  can  benefit  from 
future  software  improvements. 

3.  Software  maintenance  cost  is  low. 

4.  The  software  is  portable  to  new  platforms  and  operating  systems. 

5.  Algebraic  modeling  language's  such  as  (JAMS  allow  easy  modification 
and  addit  ion  ol  const  raints. 

We  develop  eight  tost  problems  as  described  in  Table  I.  Problem  names 
are  vague  and  data  is  hypothetical  foi  security  reasons,  but  problem  sizes 
arc-  representative  of  actual  Pakistan  Army  situations.  The  computing  time 
required  for  GAMS  to  generate  test  problems  is  given  in  Table  I.  A  large 
portion  of  the  (JAMS  time  was  taken  deriving  the  parameter  OA,//',. 

All  the  test  problems  sat isly  Condition  3  of  Section  3  for  the  initial  time 
period.  All  but  INFANTRY  I  also  satisfy  Conditions  I  and  2.  'The  test 
problems  are  all  feasible.  Tv.xpcrimenta.tion  with  other  stalling  conditions 
often  result. ed  in  infeasibility. 

The  test  problems  of  Table  1  were  solved  for  integer  solutions  with  a 
10%  optimality  tolerance  (?.c,  termination  occurs  when  the  first,  solution 
guaranteed  within  10%  of  optimal  is  obtained),  using  both  the  XA  and  OSL 
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Problem  Typo 

Problem  Size 

(JAMS  ( Jenernl  ion 
Time  (seconds) 

INFANTRY  1 

0  years 

2, 5G2  const  mints 

530 

87  unit  s 

12.  2  1 1  binary  variables 

30  locations 

74,571  nonzeros 

INFANTRY  2 

6  years 

1 ,  599  constraints 

23G 

72  units 

5,  710  binary  variables 

24  locations 

35,  037  nonzeros 

ARTILLERY 

G  years 

1,088  constraints 

151 

54  units 

3,  140  binary  variables 

21  locations 

18, 530  nonzeros 

7  years 

1 ,  5G3  constrain!  s 

200 

54  units 

5.2G3  binary  variables 

21   locations 

31 , 9G7  nonzeros 

8  years 

2, 135  constraints 

258 

54  units 

7,828  binary  variables 

21  locations 

49, 393  nonzeros 

ENGINEERS 

G  years 

913  constraints 

89 

36  units 

2,  193  binary  variables 

19  locations 

14, 829  nonzeros 

7  years 

1, 253  constraints 

116 

3G  units 

3,  421   binary  variables 

1  9  locations 

23,  37!)  nonzeros 

8  years 

J ,  684  constraints 

147 

3G  units 

4,949  binary  variables 

19  locations 

35,  130  nonzeros 

Table  1:    Test  problem  description  and  model  generation  lime  on  a   {86/33 
personal  computer. 


Problem  Type 

Solution  Time 

Herat  ions 

Horizon  (Year) 

seconds  (solver) 

INFANTRY  1 

t  =  6 

3,575  (XA) 

37,919  (XA) 

2,338  (OSL) 

7,002  (OSL) 

INFANTRY  2 

t  =  6 

1,508  (XA) 

21,047  (XA) 

•152  (OSL) 

3.  178  (OSL) 

ARTILLERY 

/  =  6 

95  (XA) 

2,519  (XA) 

173  (OSL) 

2,  190  (OSL) 

/  =  7 

079  (XA) 

14,909  (XA) 

211  (OSL) 

1,998  (OSL) 

/  =  8 

cc  (XA) 

co  (XA) 

155  (OSL) 

3,012  (OSL) 

ENGINEERS 

t  =  6 

105  (XA) 

4,735  (X A) 

105  (OSL) 

1,394  (OSL) 

/  =  7 

2,855  (XA) 

51,709  (XA) 

311  (OSL) 

3,381  (OSL) 

/  =  8 

8,950  (XA) 

105,9 18  (XA) 

1,027  (OSL) 

1  1,212  (OSL) 

Table  2:    Trsi  Problem  Solution   Times 


15 


solvers.  All  tuning  parameters  were  left  at  their  default  values.  Tabic  2 
demonstrates  these  results  and  highlights  OSL's  superior  performance  in  all 
hut  one  case.  XA  was  unable  to  solve  the  linear  programming  relaxation  of 
the  8-year  Artillery  problem  apparently  due  to  cycling. 

For  completeness  of  the  computational  testing,  all  problems  were  also 
solved  to  optimality  (i.e..  with  the  optimality  tolerance  set  to  zero).  The 
integrality  gap  (difference  between  the  optimal  linear  and  integer  program- 
ming objective  function  values)  is  zero  in  five  of  the  eight  problems  tested 
and  less  than  3%  in  the  other  cases.  For  the  five  examples  with  no  integrality 
gap,  solving  with  zero  tolerance  took  about  the  same  time  as  solving  with  a 
J 0%  tolerance;  but  the  other  three  examples  took  significantly  longer.  Since 
the  integrality  gap  was  so  small,  the  added  computational  time  for  the  zero 
tolerance  yielded  no  improvement  in  solution  quality.  Though  this  behavior 
cannot  be  guaranteed,  we  recommend  using  the  10%  optimality  tolerance 
(the  GAMS  default)  for  future  instances  of  this  problem. 

5      Conclusion 

Computational  experience  with  the  model  demonstrates  that  optimal  unit 
rotations  can  be  developed  for  6-8  year  horizons  on  a  personal  computer  in 
less  than  an  hour,  using  the  (JAMS  modeling  language  and  the  OSL  solver. 
The  Pakistan  Army  requires  five-year  rotation  schedules  for  planning  pur- 
poses, so  the  6-8  year  schedules  we  have  obtained  are  more  than  adequate 
in  scope.  The  computation  times  are  considered  acceptable  and  represent  a 
significant  improvement  ovoi  the  hundreds  of  man-hours  currently  used  to 
solve  the  problem. 

This  work  offers  some  general  lessons  for  practical  application  of  opti- 
mization modeling. 

•   Use  general   purpose  software  when   possible.     Compared   to  special- 
purpose  algorithms,  the  costs  of  development  and  long  term  mainte- 


nance  are  much  lower.    Also,  general-purpose  solvers  can  much  more 
readily  adapt  to  changes  in  I  lie  problem. 

•  Using  an  algebraic  modeling  language  to  generate  (he  model  facili- 
tates rapid  assessment  ol  computational  tracfability  and  allows  several 
solvers  to  be  tested  competitively  with  minimal  effort.  In  our  expe- 
rience, integer  programming  applications  air  too  varied  for  any  one 
solver  to  always  outperform  the  of  hois.  (E.g.,  though  XA  lost  to  OSb 
in  this  case,  it,  has  been  our  preferred  solver  in  other  applications.) 

•  When  computational  effort  prior  to  optimization  is  devoted  to  the  elim- 
ination of  unnecessary  variables,  this  refinement  can  often  make  large- 
scale  instances  of  real-world  problems  tractable. 

•  Mathematical  analysis  leading  to  additional  constraints  beyond  those 
specified  by  the  original  problem  statement,  such  as  the  constraints  on 
the  ending  conditions,  can  lead  to  better  solutions;  and,  in  the  case  of 
integer  programming,  carefully  chosen  extra  constraints  often  make  the 
model  easier  to  soke. 
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Appendix 

Theorem  If  Conditions   I,  2  and  3  of  Section  3  air  satisfied  for  six  years, 
thej'  are  sufficient  to  guarantee  a  future  feasible  rotation  schedule. 
PROOF: 

Condition  3  stipulates  that  half  the  units  in  PAs  inusl  have  previously 
served  in  II  As  and  half  in  SI  I  As.  This  status  is  maintained  on  a  yearly  basis 
provided  that  whenever  a  unit  moves  from  a  II A  to  a  PA.  anotliei  unit  moves 
from  a  SUA  to  a  PA.  The  units  al  HAs  and  SlIAs  can  therefore  be  separated 
into  pairings  which  satisfy  one  of  the  following  three  rases. 

•  CASE  1.  A  unit  has  been  at  a  II A  for  J  year  and  anotliei  unit  has  been 
at  a  SUA  for  1  year. 

•  CASK  2.    A  unit   has  been  at   a   HA  for  2  years  and  another  unit   has 
been  at  a  SUA  for  2  years. 

•  CASE  3.    A  unit    has  been  al   a   HA  for  3  years  and  another  unit   has 
been  at  SUA  for  3  years. 

It  is  shown  for  each  case  that  it  is  possible  for  the  units  at  the  HA  and 
SHA  to  feasibly  rotate  in  the  same  year  with  a  unit  currently  at  a  PA  location. 
Any  six  units  satisfying  Conditions  I,  2  and  3  for  the  last  6  years  and  feasibly 
rotating  can  be  used.  Without  loss  of  generality,  the  six  units  are  numbered 
according  to  Figure  2  (i.e..  unit  3  is  at  1 1 A  ).  For  clarify,  the  unit  idenfifiei  s  are 
bold  faced.  Also,  as  in  Figure  2,  year  0  is  considered  the  first  year  available 
to  change  a  unit's  location. 
CASE  I 

Unit  3  has  stayed  1  yeai  ;il  ;i  HA.  This  implies  unit  3  replaced  unit  f>.  I 
year  ago  and  unit  0  replaced  unit    1  either  2,  3  or  I  years  ago. 

Unit  4  has  stayed  1  year  at  a  SUA.  This  implies  unit  4  replaced  unit  5, 
1  year  ago  and  unit  5  replaced  unit  2  cither  3,  1  or  5  years  ago. 


Condition  3  ensures  that  each  year  a  move  from  a  IIA  to  a  PA  occurs. 
a  move  from  a  SUA  to  a  PA  also  occurs.  Therefore,  we  need  only  consider 
unit  6  (5)  replacing  unit  1  (2)  cither  3  or  4  years  ago. 

These  conditions  provide  the  following  rotation  eligibilities: 
Unit  4  is  eligible  to  move  in  years  1 ,  2  or  3  and  Unit   1  is  eligible  to  replace 
unit  4  in  years  2,  3  or  1  if  unit  1  was  replaced  3  years  ago,  or  in  years  1,  2 
or  3  if  unit  1  was  replaced  I  years  ago. 

Unit  3  is  eligible  to  move  in  years  0,  1  or  2  arid  Unit  2  is  eligible  to  replace 
unit  3  in  years  2,  3  or  A  if  unit  2  was  replaced  3  years  ago,  or  in  years  1,  2 
or  3  if  unit  2  was  replaced  4  years  ago. 

Therefore,  it  is  feasible  for  both  unit  3  and  unit  4  to  rotate  in  year  2. 
CASE  2 

Unit  3  has  stayed  2  years  at  a  11  A.  This  implies  unit  3  replaced  unii  G.  2 
years  ago  and  unit  6  replaced  unit   1  either  3,  4  or  5  years  ago. 

Unit  4  has  stayed  2  years  at  a  SUA.  This  implies  unit  4  replaced  unit  5. 
2  years  ago  and  unit  5  replaced  unit  2  either  4,  5  or  G  years  ago. 

Condition  3  ensures  that  each  year  a  move  from  a  11 A  to  a  PA  occurs, 
a  move  from  a  SUA  to  a  PA  also  occurs.  Therefore,  we  need  only  consider 
unit  6  (5)  replacing  unit.  1  (2)  either  A  or  -r»  years  ago. 

These  conditions  provide  the  following  rotation  eligibilities: 
Unit  4  is  eligible  to  move  in  years  0,  1  or  2  and  Unit  1  is  eligible  to  replace 
unit  4  in  years  1,  2  or  3  if  unit  1   was  replaced  A  years  ago,  or  in  years  C),  1 
or  2  if  unit  1  was  replaced  5  years  ago. 

Unit  3  is  eligible  to  move  in  years  -1,  0  or  1  and  Unit  2  is  eligible  to  replace- 
unit  3  in  years  ] ,  2  or  3  if  unit  2  was  replaced  A  years  ago,  or  in  years  0.  1 
or  2  if  unit  2  was  replaced  5  years  ago. 

Therefore,  it  is  feasible  for  both  unit  3  and  unit,  4  to  rotate  in  year  I. 
CASE  3 

Unit  3  has  stayed  3  years  at  a  11  A.  This  implies  unit,  3  replaced  unit  6,  3 
years  ago  and  unit,  6  replaced  unit  1  either  1,  5  or  G  years  ago. 

Unit  4  has  stayed  3  year  at  a  SUA.  litis  implies  unit  4  replaced  unit  5, 
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3  years  ago  and  unit  5  replaced  unit  2  either  5,  6  or  7  years  ago. 

Condition  3  ensures  thai  each  year  a  move  from  a  II A  to  a  PA  ore  his. 
a  move  from  a  SUA  to  a  PA  also  occurs.  Therefore,  we  need  only  considei 
unit.  6  (5)  replacing  unit   1  (2)  either  5  or  0  years  ago. 

These  conditions  provide  the  following  rotation  eligibilities: 
Unit  4  is  eligible  to  move  in  years  -J,  0  or  I  and  Unit  1  is  eligible  to  replace 
unit  4  in  years  0,  1  or  2  if  unit   1  was  replaced  -r>  years  ago,  or  in  years    1.0 
or  1  if  unit  1  was  replaced  (>  years  ago. 

Unit  3  is  eligible  to  move  in  years  -2,-1  or  0  and  Unit  2  is  eligible  to  replace 
unit  3  in  years  0,  1  or  2  if  unit  2  was  replaced  -r)  years  ago.  or  in  years  -1,0 
or  1  if  unit   2  was  replaced  G  years  ago. 

Therefore,  it  is  feasible  for  both  unit  3  and  unit  4  to  rotate  in  year  0.  □ 
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